AD-A179  027 


rtintnilT  Polytechnic  Institute 
School  of  Management 


DTIO 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DTIC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


GOODNESS  OF  FIT  BASED  ON 
INTEGRATED  SQUARED  ERRORS  IN 
CHARACTERISTIC  FUNCTIONS  OR  DENSITIES 


by 

A.  S.  Paulson, 
Ho-Ling  Hwang,  and 
J.  L.  Bryant 


Accession  lor 

- rr - 

NT  IS  GRAJtl 

T? 

DTIC  TAB 

n 

Un/’.nncuncsd 

□ 

Ju-tif i cntion — 

- - 

By  . - - - 

Di-tr-"  ;  '•'/ 


School  of  Management 
Rensselaer  Polytechnic  Institute 
Troy,  NY  12180-3590 
(518)266-6586 


The  view,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
authors)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 


Note:  This  paper  is  not  to  be  quoted  without  the  permission  of  the  authors). 

For  a  list  of  Working  Papers  available  from  the  School  of  Management, 
please  contact  Sheila  Chao,  School  of  Management. 


flooaw  at  Pit  Band  on  integrated  squared  Errors 
In  Characteristic  Function*  or  Densities 


by 

A.S.  Paulson* 

Bo -Ling  Hwang** 
Rensselaer  Polytechnic  institute 

and 

J.L.  Bryant 
University  at  Cincinnati 


*  Research  sponsored  In  part  by  the  U.  3 .  Army  Research  Office  under  contract 
DAA  029-81  -JC-0110. 

*•  Now  at  Oak  Ridge  National  Laboratory,  oak  Ridge,  Tennessee. 


/ 


TTtis  paper  is  concerned  aalnljr  with  tests  of  gooSsss  of  fit  thet  s  rsndosi  swpto 

*it  *tf**  **  is  from  a  cosplstsl; y  specified  distribution  function.  The  test  Is 

/ 

/ 

bawd  on  the  Integral  of  tha  watghted  aquarad  modulus  at  tha  dtfteranca  bstwssn 
aanpla  and  population  characteristic  Auctions.  This  integral  axpraaalon  ia  squivalant 
to  tha  integral  ad  tha  aquarad  diffaranca  batwaan  a  danalty  and  1M  Paraan  karnal 
MtUtete.  Tha  asymptotic  null  dlatrlbutien  ad  tha  statistic  ia  that  at  an  iidlidta 
"sighted  aa  at  Mutually  indapandant  chi  aquarad  variates.  An  approximation  to  tha 


null  distribution  ia  glvan  and  appliad  to  giwa  tha 
ad  fit  ad  p-varlate  normality.  Tha  teat  ad  Oft  ia 


points  ad  a 
mild  ragulartty 


conditions. 
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1.  introduction 


The  characteristic  function  corresponding  to  an  arbitrary  distribution  function 
P(x)  la  dafinad  by 

«..) .  r «m  lux)dP(x),  (1.1) 

where  i1  ■  -1  and  u  la  a  real  number.  If  xif  x, . .  la  a  randoai  aaagtla  froai 

P(x),  tha  sample  information  la  completely  capturad  by  tha  sample  charactarlatlc 
function 

ft 

♦n(u)  ■  n"1  £  oxpduxj)  .  (1.1) 

j«i 

tfa  ara  prlnarlly  concamad  with  tasting  tha  goodnaaa  of  fit  hypothesis 

■  i  P(x)  a  p0(x)  or  equivalently  At  a(u)  »  *s(u)  (1.3) 

against  tha  ganaral  altamatlva 

■  «  P(x)  a  P0(*>  or  equivalently  At  a(u)  *  *0(u)  . 

■ara  F0(x)  la  a  completely  apadflad  distribution  function.  Moat  of  this  papar  la 
davotad  to  tha  univariate  case.  However,  tha  structure  of  tha  analysis  and  tha 
results  claarly  Indicate  that  many  of  our  results  nay  be  extended  to  tha  p- variate 
case  with  only  minor  changes.  Discussion  of  tha  p- variate  case  is  confined  to  the 
development  of  a  test  for  p- variate  normality.  Tablee  of  percentage  points  of  tha 
statistic  and  an  approximation  to  tha  percentage  points  of  tha  asymptotic  distribution 
of  tha  statistic  ara  also  given.  Tha  distribution  of  tha  test  statistic  is  shown  to 
be  coast  asymptotic  for  vary  small  sample  sixes. 

Tha  use  of  characteristic  functions  in  tasting  hypotheses  of  fit  originated  with 
Heathcote  (1972)  who  suggested  that  the  hypothesis  of  symmetry  of  a  distribution 


V»  •»  '  .  %  •  ,  *  *  *  9  »  V»-  j.*  w 
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Z„  »  nAft  *  n  ]  l*n(u)  -  90(u)|*  ly(u)|*  du 
J  -« 


(2.1) 


»  n  r  |wn(u)  -  w0( u ) I  *  du. 

J 

sbere  o>n(u)  «  *n(u)y(u)  and  u0(u)  ■  90(u)y(u). 


(2.2) 


*•  tamtam  throughout  that  (l)  y(u)  la  tha  characteristic  function  of  mi 

absolutely  conttnuoua  distribution  function  Fy(x)  with  corresponding  daiwlty  ly(x) ,  and 
(2)  that  ly(u)  IX)  for  all  u  a  o  and  that  ly(u)  I  *  la  lntegrabla  over  ( -m,m) .  Our 
asmisfiHnne  concerning  y(u)  imply  (1)  that  all  nonsero  dlacrepandea  of  |»n(u)  • 
*0(u)  I  over  arbitrary  intervals  a<u<b,  a<b,  provide  a  positive  contribution  to  the 
integral  In  (2.1),  and  (2)  that  the  Integral  (2.1)  has  the  possibly  worm  appealing 
equivalent  fora  (Feller,  19«6,  chapter  IS,  Heathcote,  1977) 


In 


-  h0(x)>*  dx. 


(2.3) 


tew  re 


n 

h„(x)  «  n"1  £  fy(x-xj  ) 

J»i 


(2.4) 


is  an  unbiased  Par  sen  kernel  estimator  of  h0(x)  if  x,.,  xt,...,  xn  is  a  random 
sample  from  f0(x).  Since  the  convolution  of  an  arbitrary  distribution  with  an 
absolutely  continuous  distribution  is  again  absolutely  continuous  (Gukacs,  1970,  p. 

39),  h0(x)  is  the  density  associated  with  the  distribution  of  x+s,  where  x  and  s  are 
independent  random  variables,  x  has  distribution  F0(x)  and  z  has  distribution  Fy(x). 
XT  P0(x)  is  absolutely  contlnous  with  associated  density  f0(x),  then  h0(x)  * 
f0(x)*fy(x),  where  *  signifies  the  operation  of  convolution. 
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Aa  an  illustration,  vs  eonsidsr  ths  friaeworfc  for  a  tost  of  St  the  population 
distribution  et  xx,  xn  is  M(0,l).  Xn  this  cass  Se(u)  »  exp(-fcu*).  Vs 

chooss  y(u)  a  sxp(>%u*)  bscauss,  conwsnisntly,  ths  convolution  at  two  Gaussian 

variatss  is  again  Gaussian.  Vs  thus  havs 

h0(x)  ■  (2jr%)**  s*p(-%*1),  h„(x)  ■  n**(2ir)*%  C  exp( -%(x-xj  )•* 

j»x  * 

Soas  aspscts  of  ths  charactsr  of  ths  tsst  procsdurs  ars  dstsralnabls  by  snail 
seals  al Munition  experiments,  for  sxanpls  by  drawing  a  randoa  sample  xx,  xx, . . . ,  xn* 
fron  first  a  N(0,1)  population  and  aubsaquantly  fron  othsr  populations.  Plgurs 

1(a) -(c)  dapicts  hc(x)  for  this  illustration  and  typical  hn(x)  with  n«20  and  xx, 

,  Xn  a  randoa  aaapls  fron  M(0,l),  a  randoa  aaapls  froa  ths  Cauchy 

distribution  rr‘l(i+x*) *4 ,  and  a  randoa  aaapls  froa  ths  unifora  distribution  an 
(•1.14,  1.14)  rsspsctlvsly.  Ths  paraastsrs  of  ths  unifora  distribution  wsrs 

dstscainsd  so  that  its  asan,  varlancs,  and  abswnass  Matched  those  of  PQ(x) .  TMs 
figure  suggests  that  ths  tsst  for  normality  based  an  xn  will  perform  well  against 
long-tailed  alternatives  but  leas  well  against  symastric  short-tailed  alternatives,  an 
observation  verified  in  an  unpublished  Rensselaer  Polytechnic  Institute  ph.D. 
dissertation  by  Hwang  (19S4). 

Figure  l.  Comparison  of  expected  density  h0(x),  dashed,  with  h»(x),  solid, 
for  h0(x)  Gaussian  and  hn(x)  constructed  fron  (a)  Gaussian,  (b) 
Cauchy,  and  (c)  unifora  randoa  saaples. 


» 
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Since  both  l*n(u)  -  e(u)|*  and  I  y< u )  I  *  are  even  function*  of  a,  the  integral 
opn eeion  Z„  in  equivalent  to  the  Mathematically  aore  convenient 

I„  ■  n  £  ^(u)  ‘  y<u>,a  ly<u)l*  du  (2.5) 

where  the  real-valued  function*  y(u)  and  yn(u)  are  given  by 

y(u)  ■  Ree(u)  +  Za*(u),  y„(u)  »  Re*n(u)  +  Zme^u)  . 

The  function  y(u)  i*  also  a  transfers  and  yn(u)  1*  its  unbiased  sample  version. 
There  are  also  data  analytic,  graphical,  and  computational  advantage*  to  using  the 
transform  y(u)  instead  of  e(u)t  for  example  plots  of  yn(u)  and  y(u)  are  easier  to 
interpret  than  separate  plots  of  the  real  and  imaginary  parts  of  en(u)  and  e(u). 
Most  of  the  asymptotic  results  concerning  ln  are  developed  from  the  covariance  kernel 
k(u,v)  of  the  stochastic  proc**s  yn(u), 

K(u,v)  ■  n  cov(yn(u),  y„(v)) 

■  Rae( u-v)  +  Zme<u+v)  -  y(u)y(v).  (2.6) 


3.  Main  Results  Concerning  In 

Ne  now  give  the  main  results  concerning  Zn  of  (2.1)  and  (2.3).  Proof*  are 
somewhat  lengthy  and  are  not  given  here;  additional  detail*  may  be  found  ln  Bryant 
and  Paulson  (1962). 

Theorem  1.  The  goodness  of  fit  statistic  r„  »  nan  of  equation  (2.1)  and  (2.3)  has, 
when  the  null  hypothesis  H  is  correctly  specified,  the  asymptotic  distribution  whose 


characteristic  fijnctlon  c(u)  is 


where  Ax,  A,, . . .  are  the  positive  eigenvalues  of  the  Integral  operator  K  given  by 

xy(u)  ■  j*  JC(u,T)y(w)|y(v)|*  dv,  (3.2) 

where 

K(u,v)  »  R*#0(U'?)+»o(u+»)-Cl«No(uH>abo(u)JtH*«0(T)+I^0(?)]. 

This  result  characterises  the  aajnptotic  distribution  of  In  as  that  of  an  infinite 
weighted  Mia  at  Independent  chi  squared  variates,  each  having  one  degree  of  freedoa, 

i.e.,  In  has  the  asymptotic  distribution  of  £  A«x1. 

J«1 

theorea  2.  The  goodness  of  fit  test  of  the  hypothesis  Ht  e  »  #o  based  on  the 
statistic  Xn  of  (2.1)  is,  when  the  null  hypothesis  H  is  correctly  specified,  consistent. 

This  result  guarantees  that  the  power  of  the  test  approaches  unity  as  the  sasqtle 
else  increases  without  bound  for  any  alternate  P  e  P0.  Many  of  the  tests  in  the 
literature  based  on  the  use  of  the  chracterlstic  function  do  not  possess  the 
consistency  property  of  Theorea  2.  Xn  particular,  the  tests  of  Heathcote  (1972), 
Kellenaeier  (19S0),  and  Koutrouvells  (1900,  19S1),  are  all  based  on  only  a  finite 
nuaber  of  u -values  and  therefore  cannot  be  consistent. 

Theores  3.  The  Jth  cuaulant  kj  of  the  asymptotic  distribution  of  In  is 


0-1)12)-*  [A  ,  j  a  1, 

■■1  ■ 
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(3.3) 


& 


where 


The  cuaulant  expression  ( 3 .  S )  HkM  determination  of  approximations  to  the 
limiting  distribution  of  In  relatively  easy  to  obtain.  Furthermore,  this  expression  is 
easily  extended  to  dimensionality  p  >  1.  As  it  turns  out,  the  Pearson  ( 1959 )  three 
cuaulant  x*  approximation  produces  acceptable  results  In  the  upper  tall  of  the 
distribution  of  In  as  Judged  by  exact  numerical  inversion  of  (3.1)  for  a  variety  of 
special  cases.  The  approximation  determined  by  the  weighted  sum  of  two  independent 
chi  squared  variates  on  m*  and  mx  degrees  at  freedom,  axxi  +  a2x| .  Is  even  better 
but  is  a  little  more  difficult  to  compute  (Bryant  and  Paulson,  19S2). 

4.  Asymptotic  Power  of  the  Tests 

The  developments  in  this  section  closely  parallel  those  of  Durbin  and  Knott 
(1972).  Under  the  null  hypothesis  we  will  take  the  underlying  population  distribution 
Anction  P0(x)  to  be  standard  normal  with  e0(u)  ■  exp(  -yj4 )  and  will  consider  only 
simple  alternate  hypotheses  at  mean  and  variance  shift.  we  take  y(u)  = 
exp( -»§s2u* ) ,  with  a  specified  as  we  proceed .  Asymptotic  power  as  a  measure  of 
performance  of  a  test  is  not  altogether  satisfactory  since  it  is  liable  to  be  highly 
dependent  on  the  particular  null  and  alternative  distributions  chosen  for  investigation 
and  because  of  its  local  character.  Nevertheless,  the  demonstration  of  adequate 


[Wrforinct  la  terms  at  limiting  power  in  mama  given  typical  Htuation  lende  credibility 
to  the  proposed  tenting  procedure.  we  shall  omit  the  technical  detaile  at  the 
development  and  simply  pis— nt  the  results. 

Table  1  preeente  the  aeyptotic  power  at  the  goodness  at  fit  test  based  on  ln  at 
the  sx  level  at  significance  for  several  values  at  the  scaling  pan  ester  n.  Identical 
results  were  obtained  for  Ones* Hermits  quadrature  at  orders  4S  and  64,  so  it  in 
believed  that  these  are  good  approximations  to  the  asymptotic  powers  attainable 
through  une  of  Xn.  Also  Included  ln  Table  l  for  comparative  purpoeee  are  the 
asymptotic  powers  at  the  Cramer -von  Mimes  test,  Andersen* Darling  test,  and  Watson's 
0*  test,  which  have  been  computed  by  Durbin  and  Knott  (1972). 

ifote  that  the  value  at  the  sealing  factor  m  strongly  Influences  the  behavior  at  the 
test.  Large  a  yields  tests  which  are  extremely  sensitive  to  mean  deviation,  since  the 
behavior  at  the  sample  characteristic  function  near  the  origin  dominates  the  test 
statistic.  At  the  same  Urns,  the  performance  of  the  test  against  variance  shift 
deteriorates.  Conversely,  smaller  values  of  m  will,  up  to  a  point,  emphasise 
deviations  ln  variance,  if  awan  and  variance  shift  alternatives  are  at  equal  concern. 
Table  l  indicates  that  a  choice  of  m  between  .  I  and  unity  le  reasonable. 

A  comparison  at  the  ln  test  (with  m  »  l)  with  the  other  three  tests  shows  that,  if 
mean  and  variance  shifts  are  of  equal  concern,  the  Anderson -Darling  test  is  superior 
vhile  In  ranks  second.  The  asymptotic  power  of  ln  in  the  particular  situations 
addressed  here  is  therefore  comparable  to  those  of  other  commonly  accepted 
procedures.  The  attractive  power  properties  of  In  in  the  univariate  case  indicates 
that  it  Is  worth  investigating  an  extension  of  Xn  to  tests  of  p- variate  normality.  The 
tests  which  obtain  from  these  extensions  will  be  true  p-varlate  tests. 


Tabla  l 


Asyaptotlc 

Powr 

(In  S)  of  Qo 

o*»aa-of-nt  Taata 

Afainat  kmd  and  variance  Shift  in 

tha  Two-sidad  Situation 

Maan 

Shift 

varlanca  Shift 

Bast  Taat 

-  SOS 

BMt 

Taat  -  9SS 

Baat  Taat  •  SOS  Bast 

Taat 

■  ■  4/3 

47.4 

93. • 

12.7 

Sl.O 

I  ■  ■  1 

44.  • 

92.2 

17.  « 

•3.4 

n  ■  >  4/9 

41.  S 

90.3 

20.2 

•  7.3 

■  ■  2/3 

3S.4 

•a.o 

21.1 

•  7.7 

Cruar-von  Mlaas 

45.7 

92.  a 

a. s 

2S.0 

Andarson -Dari Inc 

47.2 

93.  • 

40.1 

•a. 9 

9.  Test  at  «t  for  p-variaie  (formality 


The  test  statistic  (2.1)  as  atatad  naad  not  ba  affine  invariant,  a 


drawback  in  aultivarlata  tests  at  fit  (Cox  and  Small,  1979).  Affina  Invariance  for  the 


is  however  easily  affected.  Given  a  random  sample  xt,  x. 


putative! y  from  j»p(m,D),  we  baee  a  test  of  fit  for  hp(m.D)  on 


Ip.n  ■  «  f  I9»n(u)-4<u)|*esp( -m*uru)  du,  (9.1) 


where  e(u)  ■  mxp(  -%uTu)  and 


<Mu)  ■  n'4  £  exp(  lurO'%( x.  -u ) ) . 
)««■  ■* 


(S.2) 


Here  D*^i  is  the  unique  inverse  at  the 


root  at  the  positive  definite 


covariance  matrix  o  and  u  is  a  pxi  vector  at  real  numbers.  The 


(9.1) 


has  the  dosed  form  representation 


xP,n  *  «*p  C  E  np(  -(2m* )’*<)) k) 

J  ^ 


-  2(VHS*)**9  E  9*p(  -(2(l*2ma))*1Qj  )  4-  n(l-Hiz)*%P)  (9.3) 


idiere 


Qj  -  (Xj-Ml^'^Xj-u),  Qjh  ■  <XJ-X1,)r<20)-4(XJ-XH)  . 


A  straightforward  extension  of  (3.3)  to  p -dimensions  gives  the  j-th  asymptotic 


<*  lo.n  as 


«j  ■  (j-l)i2^*1  J  Xj(  u,u)exp(  *«4uru)  du. 


(S.4) 


^  *»-»  -v  A.y.v.v 


whara  tha  Jth  itaratad  kamal  Kj (u,v)  la  daflnad  for  J  *  2  by  tha  racuralva 
{•UUcmMP 


Kj(u,v)  »  f  Kj-i<u,*»)K(w,v>OJV< -a^w)  dw 
J*P 

with  JCx(u,v)  ■  JC(u,v)  ■  <p(u-v)  -  <4»(u)<l>( -v) ,  and  whara  u,  v  and  w  ara  p*l  vactora 
oi  raal  nuabtn.  By  tadloua  Integration,  tha  flrat  thraa  cuaulanta  of  tha  aayaptotic 
dtatrtbuMon  of  Ip,n  ara  dataralnad  to  ba 

«x  -  *%»{(■*  )~fcp  -  (14«*)*%P) 

ttM  •  2lfP«a*(24aHl))_%P  -  2<(  1+0^  )*-%)"%*  ♦  (!«■*)*»)  (5.5) 

K%  •  •(a»)%P(2^t«a(34-2a*)*)‘%P  -  (3)2%pC(l*a*)(lXa^+2a4)l-*p 
♦  <3)2pC<l*a*)(H-2n*)<3+2a*)]'%p  -  C(l*a^),)‘%P). 

tha  p  dlaanalonal  varalon  of  Thaoraa  l  givaa  tha  aayaptotic  null  dlatributlon  of  Ip,n 
aa  that  of 

j><*{  • 

whara  tha  xf  v*  autually  indapandant  chi  aquarad  variatas  on  on a  dagraa  of  fraadoa 
and  tha  xp  j  ara  nomagatlve  waighta  which  aatiafy  £Ap j  *  kx. 
tha  varlata 

W  •  kx  +  (X$-v),  *  ■  »  (*••> 

*"i  Ki 

ia  conatructad  to  hava  tha  aaaa  flrat  thraa  cuaulanta  aa  tha  aayaptotic  dlatributlon  of 


X0,n»  Accordingly,  the  upper  tall  of  the  asymptotic  null  distribution  a t  xPt„  adalts  at 
Ota  ipprnil—llin  (Pearson,  1999) 

pr(  Ip t n  a  w)  »  prix6  *  v  (w«A)>  .  {9.7) 

The  approximate  null  dlatrlbutlen  of  XPf„  for  a»l  and  for  l  s  p  *  9  and  p  <  n  a 

120  has  boon  developed  by  rfulatlng  10000  independent  real  1  rati  one  of  the  atatiatlc 

XP(ft  for  each  combination  of  n  and  p.  The  cubic -apllne  — oothad  upper  10,  9,  and 

is  pointa  are  given  In  Table  2.  The  entry  under  n«a  reproeente  the  percentage  point 
computed  froa  (9.7).  Thla  table  clearly  Indlcatea  the  rapidity  at  the  convergence  of 
the  percentage  pointa  of  Xfl(P  to  thoae  of  the  aayaptutlc  dlatrlbutloa  of  Xn >p . 
Conalderable  effort  was  expended  in  an  attempt  to  find  good  valuee  of  a  aa  a  flaetton 
of  p  and  we  finally  eettled  an  the  choice  awl  aa  a  good  one,  but  It  la  not  opHaal  In 
the  aanae  of  providing  aaxlaal  power  agalnat  any  epeciflc  alternatives.  The  elhatlnn 
here  Involving  a  choice  for  a  to  provide  good  power  agalnat  unknown  aUematlvee  la 
alallar  to  that  encountered  in  density  eatlaatloni  a  variety  at  kernels  and  window 
widths  give  good  results  for  epertflc  situations  but  no  specific  choice  areas  to  be 
unlforaly  good  when  only  data,  and  not  the  parent  distribution  of  the  data,  are 
available  (Tapia  and  Thoapaon,  1979,  pp.  24-91  and  especially  pp.  40-69  and  pp. 
79-94). 

when  m  and  0  are  not  specified  and  are  estlaated  froa  the  data  by 

.  n  .  «... 

m  ■  n"1  E  *1 ,  o  ■  n'1  E  (xj-iiHxj-ar  , 

J«1  1  )■!  *  * 


define 


Z9(n(u,0)  *  n  f  IPn(u)  *  'K“)la  «»<  -a*uru)  ciu.  (5.5) 

Tha  tatUtlc  Ip,n(M.6)  hM  tha  axpUcit  capraaa ntatton  (5.3)  with  m  and  6  autoatituiad 
for  m  and  0  raapacllvaiy.  Tha  atattatlc  Ip,n(M,D)  la  afflna  invariant,  la  aa«r  to 
aamputm,  and  la  baaad  on  a  coaylata  apactOcation  of  tha  p-varlata  normal  vtdch 
axpUdtly  accounta  for  aU  davlatlena  from  p-varlata  normality  and  tta«  apadflaa  a  trua 
p-varlata  tad.  ProbaMllatle  arguaanta  atrongly  auggaat  that  tha  innrtnaaa  of  at  taat 

of  tha  ooadoalta  hypothaala  Hi  tha  parant  dlatrlbutlon  of  *4,  x, .  x*  la 

m  and  D  completely  unapadOad,  baaad  on  ZPt„(A,  D)  la  rnnalafid  for  all 
OOKmi  however,  we  hava  not  boon  abla  to  obtain  an  analytical  proof  of  tMa. 
Parcantaca  poind  for  (5.*),  l^no  ara  availahla  In  Paulaon  at  al.  (i«aa). 


AdnovftidffMinl 

Tha  Moctt  at  A.S.  Paulaon  vaa  aupportad  by  tha  O.S.  Army  Raaaarch  Office 


undar  contract  DAAC229 -ai  -K  -0110 


Bryant,  j.t.  and  Paulson, A. 3.  (1992).  Goodness  a t  fit  based  on  distance 
bet— e n  characteristic  functions  and  denaltl— .  Operations  a  a— arch  and 
Statlatlca  Research  Jtaport,  Mo.  37-92-P4,  Rensselaer  Polytechnic  institute. 

COM,  O.R.  and  small,  M.J.H.  (1979).  Tasting  multivariate  normality. 
BlomeW.  99,  pp.  293-272. 

Csorgo,  a.  (1994).  Tasting  by  tha  empirical  charactarlsttc  function*  a 
aurvay.  Proceedings  of  tha  Third  Prague  ayaposfcs  in  Asysptotlc  statistics 
(1993).  pp.  48-99.  Aastardaai  Morth-Holiand. 

csorgo,  a.  (1999).  Tasting  for  normality  in  arbitrary  H mansion.  Annals  of 
Statistics,  to  appsar. 


Durbin,  J.  and  Knott,  M.  (1972).  Conponsnta  of  Cramar-von  Mia—  Statistics.  X 
Journal  of  Ro—l  Statistics  Sodsty.  l,  34,  pp.  290-307. 

fallar,  V.  (1999).  An  introduction  to  tha  Probability  Thoory  and  Its  Applications, 
vol.  XX.  Mas  Tort i  Wiley. 

■all,  P. ,  and  walsch,  A. ■ .  (1993).  A  tost  for  normality  basad  on  tha  amptrlcal 
characteristic  function.  BlonatrUa.  70.  pp.  49S-499. 

■aathcota,  C.R.  (1972).  A  Tost  of  Ooodna—  of  Pit  for  Symaotrlc  Random  varlablas 
Austral.  Journal  at  Statistics.  £4,  pp.  172-191. 

■aathcota,  C.R.  (1977) .  Tha  Integrated  squared  error  estimation  at  parameters. 
BlometrlRa.  94,  pp.  25S-294. 

Kallarmalar,  J.  (1990).  Tha  empirical  characteristic  function  and  large  sample 
hypothesis  testing.  Journal  at  Multivariate  Analysis.  10,  pp.  79*97. 

Kautrouvells,  X.  (1990).  A  goodness -of -fit  test  at  simple  hypotheses  based  on 
the  empirical  characteristic  function.  Blometrika.  67,  pp.  239*240. 

Kautrouvells,  X.  (1991).  Test  for  normality  in  stable  law.  communications  in 
Statistics.  4,  i2,  pp.  591*599. 

tukacs,  R.  (1970).  Characteristic  functions.  Darien*  Botnar. 

Paulson,  A.S.,  Roohan,  P. ,  Hwang,  ■.,  and  filler,  M.  (1999).  A  tart  of 
goodness  at  fit  for  multivariate  normality.  Rensselaer  Polytechnic  Institute 
Technical  Report  37-99-P31. 

Pearson,  E.S.  (1999).  Note  on  an  approximation  to  the  distribution  of  non* 
central  x*.  BiometrlRa.  49,  pp.  394. 

Tapia,  R.A.  and  Thompson,  J.R.  (1979).  Monparametnc  Probability  Density 
estimation.  Baltimore*  The  Johns  Hopkins  Press. 


,  M>m-  ArJ£>.  J>mT  J-  J 


R^  M-fc  1-. 


MASTER  COPY 

_ l!MCT.ASSTPTPn _ 

gCUgTY  CL  MW  FI  C  A  now  QF  TNI*  RAGE  ?**»••  Data  bMnO 


REPORT  DOCUMENTATION  PAGE 


/\t>~/inqo2.'7 


j—^T  l  JTT  1J7TTT  T  I 


ARO-12 


«■  TITLE  tank  Subtitle) 


FOR  REPRODUCTION  PURPOSES 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


5  TYPE  OF  REPORT  *  PERIOD  COVERED 


Goodness  of  Fit  Based  on  Integrated  Squared  Errors  ,  . 

in  Characteristic  Functions  or  Densities  or  in®  PaPer 


AUTHOR*#; 

A .  S  .  Paulson 
Ho- Ling  Hwang 

.1  T  .  R  vm ran 4* 


•*  a^oihiTno  ORGANIZATION  name  ano  adoreH 

Rensselaer  Polytechnic  Institute 
Troy,  New  York  12180 


<1.  CONTROLLING  OFFICE  NAME  ANO  ADORES! 

U.  S.  Army  Research  Office 
Post  Office  Box  12211 


6.  PERFORMING  ORG.  REPORT  NUMBER 


s.  contract  or  grant  number**; 

DAAG29-81-K-0110 


>0.  program  element.  PROJECT  TASK 
AREA  •  PORK  UNIT  NUMBERS 


12.  REPORT  DATE 


>1-  number  of  pages 


MONITORINO  AGENCY  NAME  •  AOOrCSVH  BNmm  tram  Controlling  Olllca)  I  IS.  SECURITY  CLASS,  (a!  thla  report; 


Unclassified 


ISr  DECLASSIFICATION/ DOWNGRADING 
SCHEOULE 


IS  Distribution  statement  *•<  uti.  Roport; 


Approved  for  public  release;  distribution  unlimited. 


it.  DISTRIBUTION  STATEMENT  (at  the  ab  at  tact  antarab  In  Black  20,  It  til  Iterant  tram  Rapart) 


is.  supplementary  notes 

The  view,  opinions,  and/or  findings  contained  in  this  report  are 
those  of  the  author (s)  and  should  not  be  construed  as  an  official 
Department  of  the  Army  position,  policy,  or  decision,  unless  so 


■  i  i  mu  i  i  ■slwtTiir  14 


'*  REV  WORDS  (Cantlnua  an  taaaraa  alba  It  nacaaaarr  m>d  Identity  by  black  number) 

multivariate  normality,  parametric  density  estimation,  power,  generalized 
distances,  simulation,  asymptotic  distribution,  covariance  kernel,  completely 
specified  null,  characteristic  function 


- ,  mm*  HmnHty  bloc*  mwftfrj 


AMIHACr  _ _ _ _  _ _ _ 

1‘his  paper  is  concerned  mainly  with  tests  of  goodness  of  fit  that  a  random  sample 
xl>  x2»  •••»  xn  from  a  completely  specified  distribution  function.  The  test 
is  based  on  the  integral  of  the  weighted  squared  modulus  of  the  difference  between 
sample  and  population  characteristic  functions.  This  ir.terral  expression  is 
equivalent  to  the  integral  of  the  squared  iifference  between  a  density  and  its 
Parzen  kernel  estimate.  The  asymptotic  null  iistribution  cf  the  statistic  is  •'hat 
of  an  infinite  weighted  sum  of  mutually  independent  chi  squared  variates.  An 


oo un  EMIT  ON  or  I  MOV  «B  It 


OBSOLETE 


UNCLASSIFIED 


wvwww 


tecum  TV  CL*j|igic*Tigw  or  Twit  mAOtfmnm  oma 


20.  approximation  to  the  asymptotic  null  distribution  is  given  and  applied  to 
give  the  percentage  points  of  a  test  of  fit  of  p-variate  normality.  The 
test  of  fit  is  consistent  under  mild  regularity  conditions. 


UNCLASSIFIED 


JFJt- 


